Can Non-unitary Effect be Prominent in Neutrino Oscillation Measurements? 
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Subject to the neutrino experiments, the mixing matrix of ordinary neutrinos can still have small 
violation from unitarity. We introduce a quasi-unitary matrix to interpret this violation and propose 
a natural scheme to parameterize it. A quasi-unitary factor Aqf is defined to be measured in 
neutrino oscillation experiments and our numerical results show that improvement of experimental 
precision may help us figure out the secret of neutrino mixing. 
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1. Introduction The definite evidence of neutrino oscil- 
lation has confirmed that neutrinos are massive and they 
exhibit non-trivial mixing, which strengthened our be- 
lief in new physics beyond Standard Model (SM). The 
compelling experiments verifying neutrino masses and 
mixing angles are from the neutrino oscillation measure- 
ments including solar (SNO, KamLAND)[l,, ^2], atmo- 
sphere (SK)[| and reactor (CH00Z)[3 neutrino exper- 
iments. These experiments can give the squared-mass 
splitting and trigonometric function of mixing angles. Up 
to now, all the other experiments to measure the abso- 
lute values of neutrino masses can only give upper limits 
and relatively rough mixing angles. 

To explain neutrino masses, adding the neutrino 
Yukawa coupling (which gives the Dirac mass of neutrino) 
to the SM lagrangian is straightforward and causes no 
anomaly, but the real problem is the huge hierarchy be- 
tween the up component (j^g, i^^,, Vt) and the down com- 
ponent (e, /i, r) of lepton isospin doublet which is very 
different from that in quark sector. Among the mecha- 
nisms proposed to generate very light neutrino masses, 
the most popular one is the seesaw mechanism[^ 0]. In 
seesaw models heavy right-handed Majorana neutrinos 
are introduced. The Dirac masses of neutrinos at electro- 
weak scale are suppressed by Majorana mass terms to be 
ultralight Majorana neutrinos which are (primarily) left- 
handed. This mechanism embodies not only the mix- 
ing of ordinary light neutrino flavors which is similar to 
the quark mixing, but also the mixing between ordinary 
and additional heavy neutrinos which is different from 
charged fermions. 

Note current experimenters and data analyzers work in 
the framework of three-flavor neutrino, i.e., as the CKM 
matrix in quark sector, the mixing angles in parameter- 
ization of the mixing matrix extracted from the exper- 
iments just are based on a assumption that the mixing 
matrix related the neutrino mass eigenstates to the in- 
teraction ones is a 3 X 3 unitary matrix, namely PMNS 
matrixi?]. However, if there do exist more than three 
neutrino species, the matrix that transforms the ordi- 
nary neutrino mass eigenstate to the flavor eigenstate 
should not be inherently unitaryjl, Q as the additional 
heavy neutrinos can mix with the three ordinary ones 
and make the unitarity of the PMNS matrix deviated. 



The non-unitary effects of ordinary neutrino mixing 
has been extensively studied in literature, giving possible 
correction to unitarity and further with new par ameteri- 
zations for the neutrino mixing matrix [Tol. [Til. 1 1 2|. [isl. [l^ . 
While models like these emphasized on obtaining a more 
physical interpretation of neutrino mixing, there is still 
one pressing question needed to be answered: Can the 
non-unitary effects stand out within current or forthcom- 
ing measurement precision? That is to say whether the 
violation from unitarity tests in experiments will be more 
of experimental error or theoretical non-unitary effects it- 
self. The answer of the question is more fundamental in 
understanding of neutrino mixing. 

In this Letter we will start from a relatively simpler 
but reasonably parameterization for the ordinary neu- 
trino mixing matrix. A quasi-unitary matrix is proposed 
to interpret the violation from unitary. For testing the 
unitary of the neutrino mixing matrix, we will introduce 
a quasi-unitary factor Aqf , for demonstrate we also will 
give some numerical results and show in which situa- 
tion the non-unitary effects can be directly measured in 
neutrino oscillation experiments. Finally we will present 
some discussions and summary 

2. Quasi-unitary mixing and the parameterization To 
test the unitarity of neutrino mixing is important be- 
cause the unitary mixing implies only three Dirac neutri- 
nos with masses coming form the Yukawa coupling while 
the quasi-unitary mixing will imply the existence of more 
than three neutrino species and new mechanism for the 
generation of neutrino mass. Unfortunantely, up to now, 
by test of NN^ (N is neutrino mixing matrix measured 
from experiments) [9] we can not determine the error is 
whether from experiment measurements or the ordinary 
neutrino mixing matrix is non-unitary theoretically. Sub- 
ject to constraints from weak decays, the mixing matrix 
can be violated from unitarity at the order of 1% level @. 
If this mixing matrix is inherently quasi-unitary with the 
definition: 

A=(I + X)U (1) 

where X is a small matrix at subleading order (1st order), 
and U is a unitary mixing matrix: 

UU^ = I. (2) 
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If an ideal condition X + = is employed, then 

AAt=I + XX^ (3) 

thus the violation from unitarity slightly occurs at the 
second order. In this case , for three generation ordi- 
nary neutrinos, just three parameters are introduced to 
parameterize the anti-symmetric matrix X, the unitarity- 
deviation of A is at a lower order than the deviation from 
I of AAt. 

As mentioned earlier, in our letter we try to use a sim- 
ple parameterization for neutrino mixing. At first feeling, 
it seems that an anti-symmetric matrix X to interpret 
the violation is a ideal choice. However, this is not re- 
liable in many neutrino mixing models such as seesaw 
models. For example, in seesaw mechanism the mixing 
matrix for the diagonalization of left- and right-handed 
neutrinos should be an overall unitary matrix O instead 
of the 3x3 PMNS matrix for the left-handed neutrinos 
which is only part of O. To be more clear, we denote the 
mixing matrix O as 



-A-nxn Cjixn 



(4) 



where A is the matrix that transforms mass eigenstates 
of three left-handed neutrinos to the flavor eigenstates 
with n as the flavor number. The unitarity of O requires 
that: 



AA^ + CC^ = I. 



(5) 



We also denote A = (I + X)U as the quasi-unitary defi- 
nition in Eq.(IT]), where X is a small arbitrary matrix at 
subleading order. Thus we have 



AAt = I + X + X^ + XX^ = I - CCt 



(6) 



In fact, the matrix O diagonalizes the neutrino mass ma- 
trix in two steps. The first step is a block diagonalization 
which reduces the problem to a two by two problem and 
the rotation is essentially a generalized 2x2 Euler ro- 
tation. The second step is the diagonalization of light 
neutrino mass matrix by a unitary matrix Ui and that 
of the heavy one by U2 . Then the overall matrix is given 
by 



02nx2n — 



cos Ui sin Ui 
— sin U2 cos U2 



(7) 



where the sin and cos are to be interpreted as series 
expansions of a small matrix = niD /m^. From Eq.Q 



A = cos0Ui = 
C = sin0Ui = ( 



02 
2! 

03 
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4! 
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Ui, 



(8) 
(9) 



we can see that X-l-X^ is at the same order of CC^ thus 
A is violated from unitarity at the first order. Since we 



know nothing about niD and m^, in order to simplify 
the problem, we assume that 



C ^ U/r;, 

where I/77 is a small number. Then we have 



AA^ ~ I + X Xt 



(10) 



(11) 



Although the choice of the anti-symmetric X is not true 
in many models, for simplification, we still can treat X as 
a real matrix, (Strictly speaking, the matrix X should be 
hermitian in seesaw machnism, the deviation will occur 
at 1st order. Here we just use a approximate matrix to 
substitute the transformation matrix O, the diagonaliza- 
tion is not consistent with the seesaw mechanism, but it 
can simplify the parameterization for the quasi-unitary 
matrix and give us a way to test order of unitarity viola- 
tion.) thus X can be parameterized as 



X 



V 



-e X y 

-X — e z 
-y -z -e 



(12) 



where e = 1/27^. Then A is thus given by 
A2X2 = 



1 — e X 
—X 1 — e 



cos a sm t 
- sin cos ( 



(13) 



for two-flavor case and 
Asxs = 

C13C12 C13S12 Si3 

-C23S12 - S23S13C12 C23C12 - S23S13S12 S23C13 
S23S12 - C23S13C12 -S23C12 - C23S13S12 C23C13 




(14) 



for three-flavor case. Here 9 



is the rotation angle 



used to parameterize the unitary matrix U2X2 (U3X3) 
and Sij (cij ) means sin 0y (cos Oij ) . We do not consider 
the CP phase 8 in three-flavor case. In addition to the 
unitary parameterization, there are two more parameters 
for the 2x2 quasi-unitary matrix and four more for the 
3x3 one. Wc should note that 

1. Here we just treat e as a parameter which can be 
set as any small real number even zero to make A 
a quasi-unitary matrix, thus our parameterization 
is not totally based on seesaw mechanism but on a 
much simpler and more generalized way to param- 
eterize the quasi-unitary matrix; 

2. Unhke the previous works 0, [H [S [S H , in 
our parameterization AA^ is deviated from I only 
in the diagonal elements at the first order when 
e 7^ 0. This is consistent with the result of @ in 
which NNT is deviated from I apparently in the 
diagonal elements; 
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3. In seesaw interpretation, parameters e, x, y, z are 
to be too small (typically ©(IQ-^^)) to give any 
implication to the neutrino measurements, thus the 
neutrino mixing matrix A can be safely parameter- 
ized as a unitary matrix. But in general, e, x, y, z 
can be not so small, they may make unitary param- 
eterization deviated with observable errors, which 
should be theoretical errors. 

3. Implication on neutrino oscillation Now let us focus 
on the implication of such parameterization in neutrino 
oscillation experiments. Ever since the issue of neutrino 
mixing was on the table, it has always been assumed 
that the neutrino flavor eigenstates {va) (a = e, ^, t) 
were linear superposition of the neutrino mass eigenstates 
\vi) (i = 1, 2, 3) through a unitary leptonic mixing ma- 
trix: U*^ and \va) = Y^iU*i\vi). By oscillating when 
propagating in vacuum, the oscillation pro babilities can 
be expressed as (two flavors case n = 2) [l^ : 



P(lj^W) - Sapsin^ [l.27AM'^{L/E)] 



(15) 



for a ^ /3 and 



P(vl^vfi) ~ 1 - 4r„(l - T„) sin^ [l.21 /^M'^{L / E)\{1&) 

for a — f3, where E is the energy of neutrinos, L is the 
length they travel, AAf^ is the neutrino mass squared 
splitting and 



I Up 



i Up 



(17) 



Here '« up' denotes a sum over only those neutrino mass 
eigenstates that lie above AM^ or, alternatively, only 
those that lie below it. The unitarity of U guarantees 
that Sap and ATa{l — To) get exactly the same results in 
Eq.ldni) and Eq.dTO]), which yield to: 



4T„(1 -r„) = S'a/j sin^ 29. 



(18) 



This is the two-flavor case and the same result occurs for 
the three-flavor case. 

Inspired by Eq.([T8|), we define a 'Quasi- unitary factor' 
to test the unitarity of neutrino mixing: 



A 



QF 



4r,(l-T„)-5. 



For the unitary mixing Aqf — theoretically, so if we 
extract this factor from the experimental data, a non- 
zero result must have come from the experimental errors. 
But in the quasi-unitary case, a non-zero result should 
be partially from the quasi-unitarity. If e = 0, Aqf still 
equates to zero at the leading order since 



4r„(i - T„) ~ 5, 



Q/3 



sin^ 26 



2x sin 40, 



(20) 



which is similar to Eq. ((T8|) . but the situation changes 
when e 7^ 0, where we have 

4Ta(l - To) ~ sin^ 29 + 2xsinA9 - Aesin'^29 + 8esin'^9, 

(21) 

Safi ~ sin^ 29 + 2xsinA9 - Aesin^29, (22) 



therefore Aqf no longer equates to zero: 
Aqf ~ 8e sin^ 6* 0(a;^ e^ • • • 



(23) 



Thus, by comparing the neutrino oscillation data between 
P(i/q — > Va) and P{va — * J^/s), we can test the unitarity 
of neutrino mixing, one should realize the mixing angle 
9 appeared in Eq. (PT|) and Eq. ip^ has a different mean- 
ing from the ordinary one measured under the unitary 
parameterization, which is essentially an effective mix- 
ing angle, denoted as 9'^ in our subsequent discussions. 
Meanwhile, as indicated by Eq. (|^5)l . it is the diagonal 
parameter e in X that dominates Aqf while off-diagonal 
effect by x, y, z appears at subleading order thus can be 
absorbed into the effective unitary parameterization of 
the mixing matrix. But from Eq. (|20j|2^ . we know that 
x, y, z together with e can make 9'^ (maybe greatly) de- 
viated from 9 and make Aqf big enough to be measured 
in the neutrino oscillation experiments. 

Now we inspect the actual three-neutrino mixing sit- 
uation and try to give the 'quasi- unitary correction' to 
certain explicit measurements. As mentioned above, in 
order to test the unitarity we need to measure Aqf de- 
pended on a same 0^ , which needs data from different os- 
cillation patterns. In case of n = 3, however, such checks 
are very difficult, all the current accelerator, reactor, so- 
lar and atmospheric neutrino data are described within 
the framework of the 3x3 PMNS matrix. And differ- 
ent experiments measure different 9ij and Am|^- without 
considering the non-unitarity of neutrino mixing. We find 
that short distance Ve oscillation measurements may give 
us some hints. For short distance (L<5 km) it is a good 
approximation to express the Ve oscillation probabilities 
in the unitary scheme asfisj: 



1 - sin2(26li3) sm^ {/^ml^L/AE), 
sir? (26*13) sin^ (6*23) sin^ {l:\ml^L/AE) , 
sin2(26'i3) cos2(6l23) sm^ {l^ml^L / AE) . 

(24) 



^j^g") This takes the similar two-neutrino form with ^13 and 
Am§2- If the neutrino mixing is quasi-unitary, Aqf can 



be defined as: 



A 



4Te(l - T,) - Se^, - Ser 

8e8in^9f^ + Oix^, xy, xe,---), (25) 



QF 



which has the similar form to Eq. (|23p . The measure- 
ments of P{i7e Ve) have already been taken in the 
reactor neutrino experiments, however, P{ve —>■ Vfj) and 
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TABLE I: Comparison of mixing angle ranges under different 
schemes. 





Unitary scheme fl6] 


Quasi-unitary scheme 




Range 


Range 


t^l2 


31.8° - 36.4° 


24.4° - 43.1° 


nQ 
"23 


37.2° - 50.9° 


31.9° - 56.6° 


nQ 
^13 


< 10.9° 


< 18.4° 



Pii^e are difficult to measure in accelerator neu- 

trino experiments. Thus we need more precise experi- 
ments to measure Aqf for the test of unitarity. Though 
there is no data for such test, the current experiments 
data can give us some implications about the quasi- 
unitary mixing. 

4-. Testing unitarity by experiments With the expres- 
sion for the quasi-unitary factor prepared, we can present 
some numerical results and discuss the constraints on 
parameter space from current experiments and further- 
more, the probability of test the unitarity by current or 
near future neutrino experiments. As known, Solar, At- 
mospheric and reactor neutrino experiments are sensitive 
to different neutrino oscillations, which give us the best 
fit mixing angles and the errors in unitary scheme. In the 
numerical calculation, our general method was to require 
the new matrix elements Aap in Eq. (jl4p to satisfy corre- 
sponding constraints derived from the latest data[l6j. 

1) . Compare the effective mixing angles extracted in 
the unitary scheme and the quasi-unitary scheme. We 
scatter Of^ from 0° to 90°, and scan the parameters 
x, y, z randomly between —0.1 ^0.1 and e randomly be- 
tween ~ 0.1, keeping all experimentally allowed points. 
The results are listed in TABLE [T] in comparison with 
unitary scheme data derived from Ref.[l6|. 

As seen in TABLE HI with x, y, z, e 7^ and the same 
required constraints for neutrino mixing, the allowed 
ranges for neutrino mixing angles are expanded. If the 
parameters x,y,z,€ are comparable to 10~^, there will 
be highly deviated mixing angles. Though such mix- 
ing angles have no direct physical meaning (the physi- 
cal value is Aap in Eq. p^ ). we can sec from Eq. ipn]) 
that the measurement of trigonometric function has been 
changed, which means the angles we measured under uni- 
tary scheme are just effective values. 

2) . Evaluate the quasi- unitary factor in case of 
x,y,z ^ 0. Since the quasi- unitary effect in Eq. (l25l) 
can be reflected in the discrepancies of different mea- 
surements of of 2, we can also calculate the quasi- unitary 
factor from the simulation above. The numerical result 
is shown as the green points (gray x) in FIG[TJ As the 
figure indicated, under our quasi-unitary parameteriza- 
tion, nowadays experimental data cannot allow e bigger 
than 0.05 and Aqf can reach up to 2.6 x 10"^ at most, 
corresponding to deviated ^ 18.35°. 

If e and x, y, z in A are big enough, we can mea- 



sure a sizable Aqf in the quasi-unitary scheme, however, 
the error ranges in Ref.|16j] cannot be totally from the 
non-unitary effect because the experimental errors also 
take a share. Unlike the general seesaw-predicted non- 
unitary effect which is too small in comparison to the ex- 
perimental error, the quasi-unitary factor parameterized 
in this work with x,y,z ^ 0, Agp may remain sizable 
in the measurements. It is expected the precision im- 
provement of experiments will eliminate the experimen- 
tal background as much as possible, leaving a relatively 
pure constitution of Agp so that the effect may stand 
out. 




0.01 0.02 0.03 ^ 0.04 0.05 0.06 0.07 

FIG. 1: Magnitude of quasi-unitary factor versus e. The green 
points (gray x) are obtained by scattering x, y, z randomly 
between —0.1 ~ 0.1 and e randomly between ~ 0.1, corre- 
sponding to the deviation (Aqp) from unitarity theoretically. 
The blue points (dark o) are obtained by only scanning e 
{x, y, z — 0), corresponding to the deviation (Aqp) from 
unitary experimentally. 

3). Evaluate the quasi-unitary in case of x = y = z = 0. 
Since NN^ is violated from the unit matrix with defi- 
nite corrections (~ 1%) on the diagonal elements while 
only upper limits for the corrections on the off-diagonal 
elements 0, which suggests that the diagonal corrections 
dominate the unitarity violation, we also consider the 
case with a;, y, z = in our parameterization and regard 
the non-zero parameter e as being originated from the 
experimental error. This consideration is quite naive, 
though, it can help us distinguish the differences between 
theoretical unitariy violation and the experimental one. 
The simulation result (Aqp) is also displayed in FIG[T] 
by the blue points (dark o) present, for sake of comparing 
it with Aqp. 

It is clear that in this case, Agp can reach 1.0 x 10^^ 

at most, corresponding to deviated ^ 11.20°. Gen- 
erally speaking, if e in non-unitary matrix theoretically 
is the same order as the experimental error in the diag- 
onal matrix elements, Agp will be larger than Agp for 
the participation of x, y, z. This can be seen in FIG[1] 
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where Agp can be two times larger than Aqp. Mean- 
while, with the precision improvement, the tendencies of 
these two factors are to be entirely different, that Agp 
will remain non-zero while Aqp will become smaller and 
smaller and eventually vanish. 

5. Conclusion As the mixing matrix of ordinary neu- 
trinos may be unitary or quasi-unitary, we parameterize 
the mixing as a quasi-unitary matrix with four additional 
parameters e, x, y, z and define a quasi-unitary factor 
Aqf for the test of unitarity. Our numerical analysis 
gives the possible deviations in mixing angles without vi- 
olating any current experimental constraints, as well as 
the magnitude of the Agp and ASp which separately 



describes the theoretical and experimental sensibility of 
discrepancies in future ^13 measurements. The improve- 
ment of experimental precision may help us figure out the 
secret of unitarity of neutrino mixing and tell us whether 
the seesaw framework is a reasonable model for neutrino 
masses. 
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